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We propose a nonlinear extension of the 
standard tube model for semidilute solutions 
of freely-sliding semiflexible polymers. Non- 
afiine filament deformations at the entangle- 
ment scale, the renormalisation of direct in- 
teractions by thermal fluctuations, and the ge- 
ometry of large deformations are systematically 
taken into account. The stiffening response pre- 
dicted for athermal solutions of stiff rods [1] is 
found to be thermally suppressed. Instead, we 
obtain a broad linear response regime, support- 
ing the interpretation of shear stiffening at fi- 
nite frequencies in polymerised actin solutions 
[2, 3] as indicative of coupling to longitudinal 
modes. We observe a destabilizing effect of large 
strains (~ 100%), suggesting shear banding as 
a plausible explanation for the widely observed 
catastrophic collapse of in-vitro biopolymer so- 
lutions, usually attributed to network damage. 
In combination with friction-type interactions, 
our analysis provides an analytically tractable 
framework to address the nonlinear viscoplas- 
ticity of biological tissue on a molecular basis. 



1 Introduction 

The remarkable mechanical properties of animal 
cells and tissues are attributed to a dense mesh- 
work of semiflexible biopolymers known as the 
cytoskeleton. In-vitro polymerized biopolymer 



solutions provide a relatively well-controlled 
starting point for a systematic study of its ba- 
sic physical properties [4]. Particularly rele- 
vant are solutions of filamentous actin (F-actin), 
which is the universal cellular scaffolding ele- 
ment and moreover contributes to cytoskeletal 
remodelling and force generation by its tread- 
milling polymerisation [5]. While conception- 
ally simpler than crosslinked networks, pure so- 
lutions are commonly regarded as physiologi- 
cally less relevant. However, one may won- 
der about the rational basis for numerous re- 
cent "explanations" of the linear and nonlinear 
mechanics of cells and crosslinked in-vitro net- 
works as long as we do not have a sound un- 
derstanding of these properties for purely en- 
tangled solutions. Arguably, the actin bind- 
ing proteins present in the cytoskeleton of liv- 
ing cells mostly cause transient rather than ir- 
reversible crosslinking, the whole network ar- 
chitecture being constantly fiuidized by tread- 
milling and molecular motor activity. Indeed, 
numerous recent studies have established that 
the intracellular polymer network in living cells 
resembles a "soft glassy" [6,7], viscoplastic [8,9] 
material rather than a fixed elastic carcass. Its 
rheology exhibits intriguing similarities to that 
of pure actin solutions thermostated at low tem- 
peratures, which -surprisingly- show shear stiff- 
ening responses despite the absence of crosslink- 
ers [2]. Accordingly, the long-time integrity and 
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mechanical stability of the cytoskeleton seems 
to call for an explanation in terms of a merely 
physically entangled polymer solution, in which 
transient crosslinkers and other friction-like in- 
teractions affect the mechanics chiefly through 
a slowdown of the relaxation. The aim of this 
contribution is to lay the foundations for an an- 
alytical implementation of an appropriate min- 
imal model. 

The standard approach to the quasi-static me- 
chanics of entangled semiflexible polymers with 
hard core interactions employs the so-called tube 
model of semiflexible polymers — not to be con- 
fused with the tube model for flexible polymers, 
which is a dedicated scheme to deal with dy- 
namics. Its essence is a simple scaling theory 
constructed in direct analogy with the popular 
"blob model" for semidilute flexible polymer so- 
lutions [10]. In the semiflexible case, anisotropic 
blobs of length Lg and a transverse cross-section 
on the order of (P ~ Ll/£p account for the 
anisotropic fluctuations of a weakly bending rod 
of persistence length ip [11]. In the present con- 
tribution, we try to extend this successful ap- 
proach to nonlinear shear deformations. In con- 
trast to solid state physics, nonlinear deforma- 
tions are the rule rather than the exception in 
polymer networks and living cells and tissues. 
Most work done so far to understand the me- 
chanical properties of biopolymer gels has fo- 
cused on crosslinked networks composed of inde- 
pendent wormlike chains, disregarding entangle- 
ments [12-23]. For realistic fluctuating polymer 
networks, the ground state and linear mechan- 
ical properties have only very lately been ad- 
dressed on a molecular basis [24] . The nonlinear 
response was first studied in the seminal work of 
Morse extending the Doi-Edwards theory [25]. 
For filaments with lengths much longer than the 
persistence length, L ^ £p, a. softening response 
with a weak instability is obtained. The re- 
sponse to large deformations is dominated by 
"hairpins", filament segments with a large cur- 
vature in the ground state. However, the theory 
provides neither concrete results nor a physical 
picture at the mesh size scale in the biologi- 
cally relevant stiff-rod regime, L < £p. Here, 
we aim to account for the full nonlinear elas- 
tic response in a schematic 3D unit-cell of an 



entangled network of freely-sliding stiff rods, to 
obtain an analytically accessible mean-field de- 
scription in which the relevant physical mecha- 
nisms can be readily grasped. In particular, our 
approach allows for a more straightforward im- 
plementation of finite length effects and clearly 
elucidates the origin of the drastic mechanical 
difference between entangled polymer networks 
and enthalpic fiber scaffolds, which if often un- 
derestimated [1,15,16]. 

For the schematic representation of the com- 
plicated geometry of nonlinear deformations we 
follow the pioneering work by Doi and Kuzuu 
[1] for enthalpic rod networks. We borrow 
their simplified unit-cell approach, including the 
affinity assumption for the "background" of a 
test chain. This renders the problem analyti- 
cally tractable while preserving essential non- 
affinities caused by the network geometry on 
the molecular scale. The model by Doi and 
Kuzuu disregards the non-trivial physics of en- 
tanglement, which arises from the highly cor- 
related thermal fiuctuations. To include this 
crucial ingredient on an elementary level, we 
replace their enthalpic rods of fixed diameter 
by "tubes" (in the sense of the tube model 
of semifiexible polymers) of homogeneous but 
state-dependent width d. Thereby, we project 
the complicated many-body problem underlying 
the phenomenon of entanglement onto a coarse- 
grained pair problem for tubes with an effective 
pair interaction of the form of a Helfrich repul- 
sion [26]. 

Even on this schematic level, the interplay be- 
tween the free energy contributions due to bend- 
ing and tube-confinement turns out to have in- 
teresting consequences for the stress-strain re- 
lation, which hardly would have been antici- 
pated from qualitative arguments. Due to con- 
tact sliding, nonlinear bending does not play 
a significant role even for large deformations. 
Moreover, thermal undulations largely compen- 
sate the collective stiffening effect predicted by 
Doi and Kuzuu [1]. Instead, our unit-cell model 
exhibits a broad linear elastic regime. We con- 
clude that shear stiffening cannot arise in a so- 
lution of freely sliding biopolymers, supporting 
the idea that weakly adhesive interactions be- 
tween filaments are behind the responses ob- 
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served by Semmrich et al in pure F-actin solu- 
tions [2,3]. More interestingly, a geometric in- 
stability is discovered at shear deformations on 
the order of one. The corresponding "run-away" 
effects could be indicative of a shear-induced 
structural transition of the polymer solution. 
We speculate that this might then provide an 
explanation for the frequently reported collapse 
of in-vitro biopolymer solutions and networks 
under shear [2,3,27]. 

2 The model 
2.1 Overview 

The effective building blocks of an entangled 
semidilute solution of semiflexible polymers are 
thermal tubes [10]. Their free energy can be 
decomposed into two contributions: the con- 
finement of long wavelength undulations of the 
"enclosed" polymer, and tube bending. Any 
simplified analytical approach relies on a cru- 
cial assumption on how macroscopic strains cou- 
ple to the individual tube. Following Doi and 
Kuzuu [1] but replacing their rods by tubes, we 
envision the typical unit-cell test tube inter- 
acting via freely slipping entanglements with a 
background of other, identically behaving tubes. 
The test tube is assumed to be in a local equi- 
librium for the given, fixed configuration of the 
background (Fig. [T]). Thus the model describes 
the behaviour on timescales intermediate be- 
tween the terminal relaxation time and the typ- 
ical slipping time for the contacts. The fact that 
the slipping time may be very sensitive to de- 
tails of the molecular interactions and ambient 
conditions [2] raises interesting prospects for fu- 
ture developments of the model which we do not 
pursue further here. 

The major simplifying assumption of the unit- 
cell model is that the background of the test 
polymer is bound to deform affinely with the 
macroscopic deformation given by a deforma- 
tion tensor A. Due to the discrete nature of the 
contacts, the test filament deforms differently 
from its surroundings. In this way, nonaffinity 
at the entanglement scale is taken into account 
in an efficient and transparent approximation 



(Fig. [2]). The affine background has to maintain 
a constant volume to respect the (effective) in- 
compressibility of the solvent — in terms of the 
eigenvalues Aj of the deformation tensor A, 

Al A2 A3 = 1 . 

As a further, merely technical simplification [1] 
we assume filaments to lie along the eigenvectors 
of the deformation tensor with a periodic ar- 
rangement of contacts with a wavelength given 
by the mesh size. This simplification has re- 
cently been shown to introduce only very minor 
errors into the equilibrium properties of semi- 
flexible polymer solutions [24] and may be ex- 
pected to remain uncritical for our purpose. It 
has the great benefit of allowing for a simple 
intuitive explanation of the relevant physics in 
terms of a few paradigmatic cell configurations. 
Namely, the ensemble of unit-cells to be con- 
sidered reduces to the six independent realisa- 
tions depicted in Fig. [6] (see below). Consid- 
ering the full nonlinear expressions for bending 
and confinement for these representative reali- 
sations, we can then study large deformations. 

2.2 The Tube 

The free energy 

of the test tube is split into the contributions 
from confinement and bending, ^™ and A^'^, 
respectively. The confinement free energy of a 
weakly bending rod constrained by a tube-like 
harmonic potential of strength k/{2L^) is ksT 
per collision or "entanglement" length Lg [11]. 
The real constraining potential is clearly not 
uniform, but rather a discrete set of surround- 
ing tubes. Nevertheless, since the main contri- 
bution to the confinement energy comes from 
the constraint on wavelengths longer than the 
entanglement length, we may coarse-grain over 
the contact points and treat the tube as homo- 
geneous, with a (state-dependent) half width 

^/Wl:)) = d/2 ^ Ll/'i^'/' . 

For the remainder, we prefer to work in natu- 
ral units measuring energies in units of thermal 
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energy fc^T and lengths as multiples of the per- 
sistence length ^p. The confinement free energy 
of a tube segment of length ^ then takes the 
form 

And the bending energy of the space curve r(s) 
of length £ representing the tube backbone reads 



dsl 



dV 

ds2 



Taking tube backbones to be elastica [28], the 
bending energy is automatically at its minimum 
for given boundary conditions. 
We intend to describe networks of strongly en- 
tangled filaments for which both the total con- 
tour length and the persistence length are much 
longer than the mesh size ^. End-effects at the 
edges of the tube can thus be neglected. More- 
over, for our mean-field theory we may take a 
periodic arrangement of contacts with a typical 
distance ^. Since the contacts slip freely, the 
force of magnitude F acting at these contacts 
must be perpendicular to the test tube. The 
simple geometry naturally suggests the carte- 
sian coordinates illustrated in Fig. [31 The coor- 
dinate along the force is x, and the coordinate 
along the (average) tube axis is y. The point 
symmetric obstacles have (up to a sign) their 
center at Xo, Uo and contacts with the test tube 
at Xci Uc- Their relation to the backbone coor- 
dinates X£, yi at a "chemical" or arclength dis- 
tance i from the center is read off from Fig. [3l 

X£ = Xo^d = Xc^d/2 (1) 

yi = yo = yc, (2) 

the sign ± referring to the sign of yc- Because 
of the symmetry, it is sufficient to concentrate 
onto the positive sector yc> with ± — > +. 
The initial position of the contacts has an im- 
portant infiuence on the elastic response of the 
unit-cell. The mean separation between con- 
tacts along the tube is taken as the entangle- 
ment length, 2y^ = Lg. In the direction perpen- 
dicular to the tube we assume x° = d/2. Bor- 
rowing the equilibrium values for entanglement 
length and tube diameter from the analysis of 



the ground state by Hinsch et al [24] , our ground 
state conditions are 



X" 



yc 



0.32 

0.66 

+- 
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(3a) 
(3b) 



For typical in vitro actin solutions the aspect 
ratio has a quite large value of about 7 

as a consequence of the small ratio ^/Ip ~ 0.05. 
We remark that this choice of the ground state 
is only meant as a starting point. We will even- 
tually explore how much the inital aspect ratio 
affects the stress-strain relation. 
We now turn to the calculation of the bending 
energy for arbitrary states. Within our mean- 
field description it suffices to consider the prob- 
lem constrained to the plane. The direction of 
the force F is the natural reference for the tan- 
gent angle 6', given by 



dr 



Fcos{9) 



where F = |i^| is the magnitude of the external 
force (in our natural units). The equilibrium 
equation for an Euler-Bernoulli beam lying in a 
plane and in absence of body forces is 



(4) 

. Its first 
(5) 



-^=Fsme , 

the well-known equation of elastica [28 
integral is 



So-called infiexional elastica arise in our theory. 
Taking an infiexion point (where d6/ds = 0) as 
origin, the solution to Eq. O can be shown [28] 
to be given by 

cos(^/2) = k sin((/)) (6a) 
= am(s/F + /C; k) , (6b) 

where am( ; k) is the Jacobi amplitude function 
with modulus k, and the quarter-period JC is 
implicitly defined by 

am(/C; k) = 7i/2 . 
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The functional form of the solution is given by 
k, and can be therefore completely defined by 
the angle 9 = 9q between beam and force at the 
inflexion point at the origin, where 







= 7r/2 , cos(^o/2) = k , (7) 



holds. That is, the magnitude of the force enters 
merely scale factor for the arclength coor- 
dinate s (in our natural units, is roughly 
the arc length that would buckle under the load 
F). 

It will be advantageous to introduce cartesian 
coordinates (see Fig. [3]) satisfying 



^ = cos(^) 
as 



ds 



sm{9) 



f8a) 



^b) 



The solution can then be written [28] as 
2 



Xc 



S{(f);k) - S{tt/2; k) ) (9a) 



2k 



cos(0) , 



(9b) 



where £ is the elliptic function of the second 
kind. Setting initial conditions (6*, d6'/ds) at a 
given point s gives a unique solution for the elas- 
tica shape. However, the problem is rather one 
of fulfilling boundary conditions. As illustrated 
in Fig. m they are 



d^ 

— = at s 
ds 







9 = ^ 
2 



at 



(10a) 
(10b) 



The first boundary condition is already implicit 
in the form of the solution. 
A simple expression for the bending energy of 
an elastica segment can be written in cartesian 
coordinates. Integrating Eqj5] from s = to 
s = i gives 



A 



be 



cos(6^oj 



Xi 



(11) 



The tube changes its width and distorts its 
backbone to reach a free energy minimum for 
prescribed (affine) coordinates Xc, Uc of the con- 
tact point. (Alternatively, if the coordinates Xq, 



Uo are taken to obey the affine deformation, they 
have to be substituted for Xc, Uc in the following 
discussion.) 
is thus 



The free energy extremum equation 

(12) 



dA, 

^[xi,x, 

OXf 



<yc 







which 



gives the equilibrium coordinate 
xY[.Xc^ Uc) of the tube backbone. In prin- 
ciple, for the free energy minimization we thus 
have to remove F, i and 9o from A^'^ in Eq. (ITT!) 
in favour of Xc, ye. = Vc, and xe. This would 
require a full inversion of the elliptic functions 
in Eq. [TH Therefore, in practice, the derivative 
has to be rearranged into derivatives with 
respect to F, 

-1 



dA^" 




dA^" 




' dxe 




dxi 


ye 


~ dF 


ye 


(dF 


ye 1 



where all partial derivatives are taken for con- 
stant Hi. Accordingly, explicit expressions are 
only needed for xe{F, yi) and A^'^{F, yi). 
The remainder of the section is devoted to ob- 
taining expressions in terms of F, yi, Xc as inde- 
pendent variables. Setting s = £, inverting Eq. 
[9b] and replacing in Eq. [9a] gives an expres- 
sion for Xi as a function of k, F, y^. Inserting 
the second boundary condition, Eq. I10b| into 
Eqs. |6] gives 



1 

71 



k sinam(£/F + /C; k) . (13) 



The (variable) length 
squaring Eqs. [9b] and 



can be removed by 
and adding them into 



1 , ylF 

2 4 



(14) 



This equation can be used in all elliptic func- 
tions to replace the modulus k hj F and y£. It 
is then straightforward to find an analytic (al- 
beit cumbersome) expression for xe{F, ye). Sim- 
ilarly, Eqs. [13] and [7] can also be readily in- 
verted to obtain expressions for i and 9o. In this 
way we obtain the desired explicit expression 
for A^^{F, yi). Figure [5] shows dA^'^/dxely,, as a 
function of F. The two are equal for small defor- 
mations, Xi ^ ye. In this linear regime, sliding 
of the contact point is irrelevant since the lever 
arm barely changes, i.e., i ~ ye. However, as 
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the force approaches the order of magnitude of 
the spring constant, F ~ y^"^, the contact point 
begins to shde noticeably. As £ — >• oo, force and 
derivative of the bending energy uncouple: 



OA""" 



. 



yi 



Since the thermodynamic stress is given by the 
derivative of the energy, such nonlinear slid- 
ing effects can asymptotically lead to a pecu- 
liar state where the stress is zero although the 
filaments are under lateral tension, F ^ 0. 
A similar rearrangement is necessary for the 
derivative of the confinement energy, since it de- 
pends on i: 



dxp 



yc 



OF 




+ 



yc 



ddr^ 
dd/2 



where we have replaced by yc- By adding 
the derivatives of bending and confinement en- 
ergy, we finally obtain an explicit expression for 
dA/dxc as a function of F, y^. The free en- 
ergy minimization is done by numerically solv- 
ing Eq{T2] for F for the given, state-dependent 
coordinates yc and Xc- 

The ultimate goal is of course the stress-strain 
relation ct(7), with shear strain 7. To find an 
analytic expression for the shear stress cr, we 
first decompose the derivative as 



cr 



d7 



dA 

dXr 



eq 



dXr 



_^^dA 

yc d7 dye 



cq 



dye 
d7 



The equilibrium free energy can be written 
in terms of the equilibrium tube deflection 

c, yc) as 



X^c^iX 



*^(2^C) yc) — '^i^e\xc,yc)^ -^C) yc) 1 (15) 

a function of the macroscopic deformation. Tak- 
ing the derivative with respect to the contact 
position Xc gives 



dA 

dXr 



eq 



yc 



dA 

dxi 



dxr 



,yc dxc 



yc dXc 



cq 

y,xe=Xi,^ 



Since the energy is at its minimum for constant 
contact coordinates (Eq. [T^ . the first term on 
the rhs is equal to zero and therefore 



OA 

dXr 



eq 



OA 



yc 



dXr 



eq 

yc,xi=x^ 



'dd/2 



where the last equality holds since A}"^ does not 
depend on Xc- One similarly finds 



dye 



dye 



Xi 



dye 



XI 



We finally obtain 

dA'"" dxe 



a 




+ 



XI dye 




(16) 



We will later use this equation to gain physical 
insight into the mechanical stability of the so- 
lution. Note that derivatives with respect to yc 
for constant Xn require again a rearrangement, 
since we do not have explicit expressions with 
xp as independent variable. For this we need 
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dye 
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2.3 The Unit-Cell 

We have defined the x, y, z coordinates as fixed 
to an arbitrary set of tubes. The test tube hes 
along y and interacts with background tubes 
which lie along z and move in the x-y plane. 
For an isotropic network, the unit-cell should be 
modeled superposition of test- and back- 
ground tubes lying along all possible orienta- 
tions, integrating the triad x-y-z over the whole 
3D triad-space. However, this would compli- 
cate the problem enormously, since triads ly- 
ing oblique to the main axis of the deforma- 
tion would distort. To simplify matters we fol- 
low Doi and Kuzuu [1] and distribute test- and 
background tubes only along the main stretch 
directions (the eigenvectors of the right stretch 
tensor [29]). Since then tubes remain perpen- 
dicular to each other, the deformation can be 
incorporated in the analytical treatment as a 
translation of the contact point. The total free 
energy of the unit-cell is then a sum of 6 terms, 



A = Ai2+ ^13 + A21 + ^23 + ^ 



^31 



A 



32 
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We denote configurations by the coordinates of 
the deformation plane; for a particular tube 
triad, the free energy is given by Axy E Figure 
[6] depicts all possible arrangements. 

2.4 The deformation 

We assume the contact points to deform affinely, 
going from their initial positions Rc to final po- 
sitions Tf.. The assumption of affinity means 
that the deformation gradient A = drc/dHc, 
is spatially constant, and thus 

= A Rc 

holds for all contact points. With our assump- 
tion that the tubes lie along the stretch direc- 
tions, this statement can be written in terms of 
the stretch factors as 



Vc = ^Vc ■ 

Solving Equation [12] we obtain the equilibrium 
free energy A'^^ ClS Sb function of the deforma- 
tion. The (Cauchy) stress tensor cr can then be 
calculated as the derivative of the equilibrium 
free energy per unit volume, 

cr = detA — — ^ — A , 
aA 

where the volume of the unit-cell is taken as 



0„,0 2 



y, 



For concreteness and comparison with experi- 
ments we will restrict ourselves to volume pre- 
serving shear deformations in the plane 1-2, 
with positive stretch along 1 and contraction 
along 2. In terms of the eigenvalues of the strain 
tensor, Ai ^ oo, A2 ^ 0, and A3 = 1 when 
7 00. The shear stress is given by 

dA'^'^/Ve 



a 



(17) 



Since the ground state is assumed to be the 
same for all orientations — as expected for an 
isotropic network — the sum of all stresses must 
give a = for 7 = for the unit-cell. For shear 
in the 1-2 plane, it is readily seen that two con- 
figurations cancel each other's stress if they are 
related by an interchange of the directions 1, 2. 



2.5 Doi-Kuzuu effect 

By "Doi-Kuzuu effect" we mean the change in 
the number of contacts between filaments due 
to the deformation [1]. A parallelepiped with 
sides 2d/\x, L/Xy, L/X^, assumed to lie along 
the principal stretches, distorts affinely into a 
parallelepiped with sides 2d, L, L. By con- 
struction, the number of filaments lying along 
z which come into contact with the test rod can 
be readily shown to be [1] 



An 



c 



"i XyXz \ Xx 



2d 



- 2d 



cdL^il - A, 



Thus, as the solution is distorted at constant 
volume, the number of contacts between rods 
increases, reducing the lever arm for bending. 
Assuming stress to be of purely enthalpic nature 
leads to a pronounced stiffening response [1]. 
We remark that the material becomes stiffer 
without the rods ever leaving the linear bending 
regime, a pure network effect. 
Such a general stiffening mechanism is certainly 
worth considering and can be readily incorpo- 
rated in our theory. The inverse of the number 
of contacts per unit length, L/n, amounts to 
the distance Uc between contacts. For our pur- 
poses, however, the stiff rods have to be replaced 
by thermal tubes with a variable diameter. For 
this we recast the original argument in differen- 
tial form. Replacing cL/3 by 1/^^, we get 



he 



2_yld 

e 



SX. . 



An extra term is needed to address tube thick- 
ening in absence of shear, which increases the 
number of contacts according to 



2y| 



Sd 



The two terms influence the evolution of the dis- 
tance between contacts, yc- Within our mean- 
fleld approach it is natural to add them to the 
affine macroscopic deformation: 



*Here we deviate from the notation used in Ref. [1] 



Syc = yc SXy + 



Mi?«,_!|i,. (18) 
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3 Results 

3.1 Softening and instability 

Figure [7] shows the response of the unit-cell for 
mesh sizes ^/d-p = 0.01 and 0.1. The stress- 
strain relation shows a broad linear regime 
which softens and becomes unstable beyond 
100% shear strain. Though the theory employs 
the full nonlinear equations for bending, stiff- 
ening is not observed, suggesting that the de- 
flection of the test tube remains small through- 
out. Indeed, linearizing A^'^ in xi around zero 
has only a minor effect on the stress-strain re- 
lation, as shown in Fig. [71 A consequence of 
the small deflection is that changing the mesh 
size ^ essentially amounts to a scaling of the 
stress axis, as can be seen in Fig. [71 For lin- 
earized bending the equilibrium deflection scales 
like x'p ~ which renders both confinement 
and (linearized) bending energy independent of 
the mesh size (both are ~ ksT). Therefore the 
stress goes like 
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which boils down to the well known prediction 
of the tube model for the concentration depen- 
dence of the shear modulus, Gq ~ c^^^. Devi- 
ations from this simple mesh-size scaling are a 
measure of nonlinear bending. 
A more important parameter is the initial as- 
pect ratio, Vc/x^', as shown in Fig. [HI it decides 
on the stability of the response. A unit-cell with 
an aspect ratio larger than ~ 10 is already un- 
stable in the ground state. Smaller aspect ratios 
push the instability away, but without a signif- 
icant increase in the modulus. 
To better understand the mechanical response 
of the unit-cell in terms of the microscopic de- 
formations, it is helpful to consider the individ- 
ual contributions of each xy orientation (Fig. 
[H]). In general, none of the individual configura- 
tions displays stiffening behaviour. To a varying 
extent the instability shows up in all contribu- 
tions. It appears at its strongest in the 2-1- 
term, corresponding to a test tube along the 
stretch direction and background tubes along 
the third direction, moving in the plane 1-2 (see 



Fig. [H]). Since the transverse stiffness of a beam 
of length L decreases proportional to with 
increasing length, an increase in the mean dis- 
tance between adjacent contacts i/c weakens the 
resistance of the test tube, so that its backbone 
yields to the pressure exerted by the surround- 
ing tubes, which try to expand in order to min- 
imize their confinement free energy. As the de- 
fiection of the test tube increases, the negative 
contribution to the stress from 

dye 



will eventually dominate Eq. [161 The instabil- 
ity is an inescapable consequence of our Ansatz 
for the background deformation and the freely 
sliding contacts. Figure [TOl tries to convey an 
intuitive feeling for it. 

One may wonder whether this instability is an 
unavoidable feature of stiff polymer solutions 
or whether higher order correlations neglected 
in our simple unit-cell approach will eventually 
stop a real solution from yielding. After all, 
a living cell should be able to withstand large 
stretch without collapsing. In the remainder of 
this section we will digress on two mechanisms 
which provide stability against large deforma- 
tions: the Doi-Kuzuu effect [1] and hairpins [25]. 

3.2 Thermal Doi-Kuzuu effect 

The common trend in the results shown so far 
is a softening response. This seems at odds 
with the experimentally observed stiffening re- 
sponses in F-actin solutions [2,3]. The question 
arises whether collective effects may be behind 
them, as in the Doi-Kuzuu effect, which pre- 
dicts a power-law stiffening response for a solu- 
tion of rigid rods. As shown in Fig. [Til replac- 
ing the rods by tubes with a variable diameter 
has a dramatic consequence: the stiffening re- 
sponse becomes a softening one. The tug-of-war 
between Doi-Kuzuu reinforcement via increase 
in entanglement density and the thermal slid- 
ing instability leads to a neutral result where 
both largely compensate each other. Therefore, 
though unable to provide stiffening, the Doi- 
Kuzuu effect renders the unit-cell stable. 
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3.3 Hairpins are stable 

Here we point out a simple mechanism to sta- 
bilize the network even on the level of the sin- 
gle unit-cell. Namely, as depicted in Fig. O if 
the contact Xc lies "on the other side", Xc > 0, 
such that the test tube is more strongly bent 
in the ground state. This boils down to a mi- 
croscopic realization of the "hairpins" first dis- 
cussed by Morse [25]. Though at first sight in- 
nocuous, the change in topology has dramatic 
consequences for the stress-strain relation. As 
shown in Fig. [13], the response becomes a weak 
stiffening one and the sliding instability van- 
ishes. This weak stiffening is not a consequence 
of nonlinear bending, as can be seen by the re- 
sponse of the linearized theory. Stiffening arises 
as the distance between contacts and hence the 
lever arm i/c decreases, an effect which domi- 
nates for the configuration 1-2 (see Fig. [131 
inset). The nonlinear theory actually tames the 
response into softening at very large shear. This 
is a consequence of entering a purely sliding 
regime as in Fig. [5l Thus, though the linearized 
theory suffices for the typical unit-cell, it grossly 
overestimates stiffening in the hairpin configu- 
ration. 

4 Discussion 

In the present work we have developed a sim- 
ple, analytically tractable mean-field descrip- 
tion of solutions of stiff biopolymers, describ- 
ing an adiabatic evolution at timescales longer 
than typical times for crosslink slippage. The 
theory makes use of the full nonlinear bend- 
ing equations and is hence suitable to address 
large deformations. Surprisingly, in spite of 
all stiffening mechanisms considered, the the- 
ory invariably predicts a broad linear regime 
which softens and eventually becomes unstable 
at large shear strain. Nonlinear bending, which 
leads to strong stiffening responses in the fixed- 
contact scenario [30], becomes largely irrelevant 
when contacts slide. The hairpin configuration, 
though stable, gives a very weak stiffening re- 
sponse which does not resemble the experimen- 
tal data, where the linear regime does not ex- 
tend beyond 5% strain [2,3]. Finally, our ver- 



sion for semiflexible polymers of the geometric 
Doi-Kuzuu stiffening mechanism [1] is largely 
cancelled out by thermal fluctuations. Though 
it does provide stability, the response remains a 
softening one. We conclude that stiffening in 
biopolymer solutions necessarily requires cou- 
pling to the longitudinal stretching and com- 
pression modes of the individual filaments. Our 
conclusion supports the previous interpretation 
of stiffening in terms of the glassy wormlike 
chain [2]. Understood as a consequence of fric- 
tion, the fact that stiffening requires fast rates 
and low temperatures and leads to large dissi- 
pation makes indeed sense. 
Having gained this insight, we can outline a full 
theory of nonlinear viscoelasticity of biopolymer 
gels. At fast rates, effectively the solution be- 
haves like a network; here one may expect theo- 
ries developed for crosslinked networks to hold. 
Importantly, the elastic strain will be limited to 
values below ~ 5% strain, since at larger strains 
crosslinked networks stiffen dramatically. The 
large forces attained by the stiffening filaments 
will break bonds faster and render the deforma- 
tion inelastic. Thus, for finite rates we expect 
the nonlinear response to go in hand with huge 
structural dissipation. As the rate is lowered to 
timescales intermediate between crosslink slip- 
page and terminal relaxation, the present theory 
should hold, predicting an almost reversible re- 
sponse, with a broad linear regime. This picture 
is indeed corroborated by experiments. Both 
for living cells [8,9,31,32] and in-vitro F-actin 
solutions [2, 3], stress-strain relations become 
increasingly linear over large deformations as 
the strain rate is lowered, and the reversibil- 
ity of the response increases as dissipation goes 
down. Thus, there seems to be hope of a uni- 
fied description of this plethora of mechanical 
responses in terms of networks of semiflexible 
filaments with frictional interactions. 
As the deformation spreads the contacts along 
the main stretch axis, their lever arm increases 
and the test tube no longer can sustain the ther- 
mal fluctuations. This leads to a sliding insta- 
bility and presumably to shear banding at a 
critical strain of 100%. This picture provides 
a fresh, novel look at the mysterious "collapse" 
of biopolymer networks and solutions. Pure F- 
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actin solutions undergo irreversible weakening 
at ~ 150% [2,3]. The experimental fact that 
this collapse takes place at a well defined strain, 
largely independent of rate and stress, makes 
the case for a purely geometrical instability, un- 
related to breakage of filaments, possibly resem- 
bling collective rearrangements in oil tubes [33]. 
Though the full nonlinear theory including the 
coupling to the longitudinal polymer modes is 
still under development [34] , the reasoning out- 
lined above supports the sliding stability 
qualitative explanation for the collapse. Indeed, 
since the sliding instability takes place at strains 
~ 100% in our theory, and the backbone strain 
cannot be much larger than 5%, one would ex- 
pect the critical collapse strain to be largely un- 
affected by the degree of stiffening, i.e., by the 
strain rate - as experimentally observed. 

Hairpin-entanglements are interesting, as they 
are stable against large deformations. In fact, 
these hairpins are always present with a certain 
small probability on the scale of the unit-cell, 
and with higher probability on larger scales. 
Closely related to the hairpins described by 
Morse in long filaments [25], they have a dras- 
tic effect onto the mechanical response, because 
their deformation and free energy invariably 
increase upon affine background deformations. 
By addressing hairpins and normal entangle- 
ments separately, their opposite responses at 
large strains stand out. The instabihty of the 
unit-cell against shear may be turned into a 
weak stiffening response upon changing the en- 
tanglement topology to a hairpin. One may 
therefore wonder whether hairpins are behind 
the fact that the conspicuous in-vitro collapse 
is not observed when straining living cells by 
over 100% [8,9]. Since large deformations of 
this order are to be expected under physiolog- 
ical conditions (as cells crawl and divide), it 
seems desirable to avoid a collapse of the main 
load-bearing structure in the cell. Though typ- 
ical actin filament sizes are probably too short 
to provide strongly bent tubes, one may cre- 
ate hairpins in alternative ways. To achieve 
the right "inside-out" topology, it would suffice 
to crosslink filaments at a large angle and let 
them polymerize, thereby effectively generating 
a large bent in the ground state. Such a strategy 



would remind of the Arp 2/3 complex, an ubiq- 
uitous component of the actin cortex. As a more 
generic possible mechanism one can imagine a 
pre-bending of filaments in the ground state en- 
forced by mutual adhesive interactions or sim- 
ply by their embedding into a disordered envi- 
ronment. This would provide an intimate link 
between the question of mechanical stability of 
living cells and tissues and a fundamental poly- 
mer physics problem. 

We are very grateful to Andreas R. Bausch 
for his generous support. Insightful discussions 
with C. Semmrich and J. Glaser are also ac- 
knowledged. 
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Fig. 1 The test tube configuration is a com- 
promise between bending and con- 
finement free energy for fixed posi- 
tions of tlie contacts witli tlie sur- 
rounding tubes. 

Fig. 2 Tlie background deformation is as- 
sumed to be afiine and volume- 
preserving. The tube adapts adi- 
abatically to its local free energy 
minimum. Note that the tube de- 
forms non-affinely although the con- 
tact points obey the macroscopic 
(affine) deformation. The number 
of contacts may also change as the 
background is distorted [1]. 

Fig. 3 A tube interacting with the sur- 
rounding tubes. The tube width is 
d and its deflection is xi. Taking ad- 
vantage of the simple geometry, the 
cartesian coordinate system is cho- 
sen with X along the force F and y 
perpendicular to the force along the 
(average) tube axis. Two alternative 
versions of the central assertion are 
considered, namely that (1) the con- 
tact point (black dot) coordinates Xc, 
He or (2) the obstacle centers Xo, Uo 
vary affinely with the macroscopic 
strain. 

Fig. 4 Boundary conditions for the tube 
backbone. At the origin, s = 0, 
symmetry requires zero curvature, 
dO/ds = 0. Since the entanglement 
at s = £ slips freely, the force F must 
be perpendicular to the tangent; i.e., 
ei = 71/2. 



Fig. 5 



Derivative dA^'^ / dxe\y^ as a function 
of force F at the contact point (in 
units kbTip/y'^). Only in the lin- 
ear regime they can be identified. 
For large deflections X£ ^ ye, the 
force saturates at 2kBTip/yj and the 
bending energy no longer grows with 
the backbone coordinate Xf. 



Fig. 6 Unit-cell for shear in the 1-2 plane. 

The main stretch directions 1-2-3 



rotate with the current configura- 
tion. The unit-cell is modeled as 
a superposition of test tubes lying 
along the main stretches. Each con- 
figuration is named according to the 
x-y coordinates. The arrows indi- 
cate the displacement of the contacts 
upon an increase of the affine strain 
7- 

Fig. 7 Normalized shear stress-shear strain 
relation for mesh sizes ^/ip = 0.01 
(sohd lines) and 0.1 (dashed lines). 
The mesh size barely makes a differ- 
ence, a consequence of the small fil- 
ament defiection. The deformation 
of the background has been consid- 
ered in two ways. Gray lines: middle 
points deform affinely, Xq = XxXq. 
Black lines: contact points deform 
affinely, Xc = XxX^- Dotted line: 
linearized bending energy. For the 
linearized theory, changing the mesh 
size merely scales the stress axis. 

Fig. 8 Aspect ratio dependence. Shear 
stress-shear strain relation for mesh 
size ^ = 0.02 ^p. The initial aspect 
ratio ?/c/^c increases along the ar- 
row, taking the values 3.5, 4.9, 7, 10, 
14, with ~ 7 corresponding to 

the tube model, Eq. [31 

Fig. 9 Shear stress versus shear strain for 
all 6 possible orientations of the test 
and background tubes, for mesh size 
^ = 0.02 £p. The curve passing 
through zero stress is the sum of all 
others, representing the response of 
the unit-cell. The instability is es- 
sentially due to the term 21. 

Fig. 10 Cartoon illustrating the instability 
by an intuitive analogy. Due to in- 
compressibility of the shaded area, 
the contact point (represented by a 
filled circle) can only move along the 
dashed line. The springs (represent- 
ing confinement free energy) have an 
infinite rest length and slide along 
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clamped beams (representing tube 
bending). Movement of the contact 
point increases (lowers) the compres- 
sion on the spring opposing (favour- 
ing) the movement. The arrows cor- 
respond to the total force acting on 
the contact point. In the ground 
state the lever arm is short and the 
beams barely bend - the system is 
linearly stable. As the network is 
sheared further the increase in lever 
arm softens the beam, whose bend- 
ing gives rise to a lateral force push- 
ing out the contact. The system be- 
comes unstable. 

Fig. 11 Doi-Kuzuu effect. Shear stress-shear 
strain relation for mesh size ^ = 
0.02 £p, with and without DK effect. 
The stiffening response reported in 
Ref. [1] is "ironed out" by the ther- 
mal fluctuations. 

Fig. 12 A hairpin. Setting a;° > 0, while 
still complying with the requirement 
kcl be I — leads to a very different 
response. Since now the tubes are 
"trapped", the entangled solution is 
stable against large deformations. 

Fig. 13 Shear stress versus shear strain for a 
hairpin configuration where Xc > 0, 
as in Fig. [121 Solid line: full nonlin- 
ear theory. Dashed line: linearized 
theory. The instability vanishes for 
the hairpin configuration. Note that 
the weak stiffening response is not 
due to nonlinear bending, since it is 
present also in the linearized theory; 
it is rather a consequence of the de- 
crease of the distance between con- 
tacts. Inset: Shear stress versus 
shear strain for all 6 possible orien- 
tations of the test and background 
tubes. 
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Figure 1: The test tube configuration is a com- 
promise between bending and confinement free 
energy for fixed positions of tlie contacts witli 
the surrounding tubes. 
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Figure 2: The background deformation is as- 
sumed to be affine and volume-preserving. The 
tube adapts adiabatically to its local free energy 
minimum. Note that the tube deforms non- 
affinely although the contact points obey the 
macroscopic (affine) deformation. The number 
of contacts may also change as the background 
is distorted [1]. 
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Figure 3: A tube interacting with the surrounding 
tubes. The tube width is d and its deflection is 
xi- Taking advantage of the simple geometry, the 
cartesian coordinate system is chosen with x along 
the force F and y perpendicular to the force along 
the (average) tube axis. Two alternative versions 
of the central assertion are considered, namely that 
(1) the contact point (black dot) coordinates Xc, He 
or (2) the obstacle centers Xq, Ho vary affinely with 
the macroscopic strain. 
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Figure 4: Boundary conditions for the tube back- 
bone. At the origin, s = 0, symmetry requires zero 
curvature, dO/ds = 0. Since the entanglement at 
s = I sHps freely, the force F must be perpendicu- 
lar to the tangent; i.e., = tt/2. 
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Figure 5: Derivative 9^ /Ox^l^^^ as a function of 
force F at the contact point (in units k},TlplyJ). 
Only in the linear regime they can be identified. 
For large deflections yi, the force saturates at 
2kBTip/yj and the bending energy no longer grows 
with the backbone coordinate X£. 
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Figure 6: Unit-cell for shear in the 1-2 plane. The 
main stretch directions 1-2-3 rotate with the cur- 
rent configuration. The unit-cell is modeled as a 
superposition of test tubes lying along the main 
stretches. Each configuration is named according 
to the x~y coordinates. The arrows indicate the 
displacement of the contacts upon an increase of 
the affine strain 7. 
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Figure 7: Normalized shear stress-shear strain 
relation for mesh sizes ^/d-p = 0.01 (solid 
lines) and 0.1 (dashed lines). The mesh size 
barely makes a difTerence, a consequence of 
the small filament deflection. The deforma- 
tion of the background has been considered in 
two ways. Gray lines: middle points deform 
affinely, Xq = XxXq. Black lines: contact points 
deform affinely, Xc = Xx^c- Dotted line: lin- 
earized bending energy. For the linearized the- 
ory, changing the mesh size merely scales the 
stress axis. 




Figure 8: Aspect ratio dependence. Shear 
stress-shear strain relation for mesh size ^ = 
0.02 ip. The initial aspect ratio increases 
along the arrow, taking the values 3.5, 4.9, 7, 
10, 14, with ~ 7 corresponding to the 

tube model, Eq. [31 
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Figure 9: Shear stress versus shear strain for 
all 6 possible orientations of the test and back- 
ground tubes, for mesh size C, = 0.02 ip. The 
curve passing through zero stress is the sum of 
all others, representing the response of the unit- 
cell. The instabihty is essentially due to the 
term 21. 



Y«l, 

linearly stable 




Figure 10: Cartoon illustrating the instability 
by an intuitive analogy. Due to incompressibil- 
ity of the shaded area, the contact point (repre- 
sented by a filled circle) can only move along the 
dashed line. The springs (representing confine- 
ment free energy) have an infinite rest length 
and slide along clamped beams (representing 
tube bending). Movement of the contact point 
increases (lowers) the compression on the spring 
opposing (favouring) the movement. The ar- 
rows correspond to the total force acting on the 
contact point. In the ground state the lever arm 
is short and the beams barely bend - the system 
is linearly stable. As the network is sheared fur- 
ther the increase in lever arm softens the beam, 
whose bending gives rise to a lateral force push- 
ing out the contact. The system becomes un- 
stable. 
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Figure 11: Doi-Kuzuu effect. Sliear stress-sliear 
strain relation for mesli size C, = 0.02 ip, witli 
and witliout DK effect. Ttie stiffening response 
reported in Ref. [1] is "ironed out" by tlie tlier- 
mal fluctuations. 



25 




Figure 12: A "hairpin". Setting x° > 0, while 
still complying with the requirement ~ 
leads to a very different response. Since now 
the tubes are "trapped" , the entangled solution 
is stable against large deformations. 
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Figure 13: Shear stress versus shear strain for a 
"hairpin" configuration where > 0, as in Fig. 
[T2l Sohd fine: full nonlinear theory. Dashed 
line: linearized theory. The instability vanishes 
for the hairpin configuration. Note that the 
weak stiffening response is not due to nonlinear 
bending, since it is present also in the linearized 
theory; it is rather a consequence of the decrease 
of the distance between contacts. Inset: Shear 
stress versus shear strain for all 6 possible ori- 
entations of the test and background tubes. 
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